Interaction Effect in the Kondo Energy of the Periodic Anderson-Hubbard Model 
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We extend the periodic Anderson model by switching on a Hubbard U for the conduction band. 
The nearly integral valent limit of the Anderson-Hubbard model is studied with the Gutzwiller 
variational method. The lattice Kondo energy shows [/-dependence both in the prefactor and the 
exponent. Switching on U reduces the Kondo scale, which can be understood to result from the 
blocking of hybridization. At half-filling, we find a Brinkman-Rice-type transition from a Kondo 
insulator to a Mott insulator. Our findings should be relevant for a number of correlated two-band 
models of recent interest. 



The usual scenario of the Kondo effect (either for a 
single impurity, or for the Kondo lattice) envisages that 
the localized moments of strongly correlated electrons are 
spin-compensated by a non-interacting electron gas. Re- 
cently, there has been a great deal of interest in the more 
difficult problem when the spin-compensating medium 
itself is a more-or-less strongly correlated electron sys- 
tem characterized by its own Hubbard U. Attention has 
been concentrated on the impurity case (!]], but there 
are also preliminary results on periodic models such as 
the Kondo-Hubbard lattice model ||. In all studies, one 
finds a strong [/-dependence of the Kondo energy. The 
subsystem described by the Hubbard model can be re- 
placed by other correlated models such as the t-J model 
H , the quantum Heisenberg model Q, or a Luttinger 
liquid §. 

The standard Kondo problem can be treated at two 
levels: one can start either from the purely fermionic An- 
derson model, or from the Kondo Hamiltonian which is 
the large- E/ effective Hamiltonian of the Anderson model. 
In the relevant parameter range the two approaches lead 
to equivalent results. Remarkably, once the Hubbard U 
of the conduction electron sea is switched on, it mat- 
ters very much whether the impurity is described by 
the Anderson or the Kondo model: one finds a reduced 
Kondo coupling for the Anderson impurity [Q , while the 
Kondo energy is enhanced if the Kondo Hamiltonian is 
used H . For an Anderson impurity in a Luttinger liquid, 
one even finds a non-monotonic [/-dependence, with the 
Kondo temperature increasing at weak U, and decreas- 
ing at strong U ||. These findings are not necessarily 
mysterious: the handwaving argument might be that in 
the Anderson case, switching on U blocks hybridization 
and therefore decreases the effective Kondo coupling. In 
contrast, in the Kondo Hamiltonian the coupling is al- 
ready given, and then switching on U helps the singlet 
formation because a more correlated electron gas offers 
more uncompensated spins. It appears that the Kondo 
impurity model with a Hubbard band is not the effec- 
tive Hamiltonian for an Anderson impurity in the same 
Hubbard band, and it is worth to explore the two prob- 
lems independently. The same holds for periodic mod- 



els where it depends on the physical nature of the sys- 
tem whether the Anderson-Hubbard lattice model, or a 
Kondo-Hubbard lattice model is a better starting point. 
We note the suggestion [Q that the curious correlated 
semiconductor FeSi can be described by an Anderson- 
Hubbard-type model in which the spin-compensating and 
spin-compensated electrons play symmetrical roles. Per- 
haps a Kondo-typc Hamiltonian gives a better starting 
point for the Nd2- a: Ce a; Cu04 system J|||. 

We consider the periodic Anderson-Hubbard model 
which describes an array of strongly correlated /-sites hy- 
bridized with a moderately strongly interacting d-band: 
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where h{ — u. fin-, etc., the k are wave vectors, and the 

JCT J JfT J J u ' 

j are site indices. The <i-bandwidth is W. In what fol- 
lows, we take the strongly asymmetric Anderson model 
with Uf — > oo and the /-level e/ < sufficiently deep- 
lying so that we are in the Kondo limit: 1 — nt <C 1 



where the /-valence is defined as nt 
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total electron density (per site, for one spin) is n = 
(V. h( + V. ht r )/2L, where L is the number of lat- 
tice sites. We will assume 1/2 < n < 1, so that there 
are enough electrons to fill at least the /-levels, and the 
d-band filling is variable up to half-filling. 

We use the Gutzwiller variational method, generaliz- 
ing a previous treatment of the periodic Anderson model 
]i~0| |. Limitation of space preventing us from going into 
details, we just mention a few salient points. The trial 
state is obtained by Gutzwiller-projecting the optimized 
hybridized-band state 
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where the mixing amplitudes itk/fk are treated as inde- 
pendent variational parameters. The Gutzwiller projec- 
tor for the d-electrons is 
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where the variational parameter rj is controlled by Ud- 
For the /-electrons the full Gutzwiller projection is taken 
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Here we consider only non-magnetic solutions corre- 
sponding to a mass-enhanced metal at n < 1, or a 
renormalized- hybridization-gap insulator at n = 1. The 
study of magnetic ordering, especially in the less-than- 
half-filled case, should prove interesting: metallic fer- 
romagnetism has been found in the Kondo-i-J model 
H , as well as in our preliminary investigation pl| of the 
Fu-Doniach model 0. We note that Tasaki's two-band 
model for which (for a certain choice of the parameters) 
the existence of ferromagnetism can be shown exactly 
Jl2| , is also of the form (]]]) but with a dispersive /-band 
and a k-dependent hybridization. 

The lengthy details of the optimization procedure, 
along with results for the behaviour away from the nearly 
integral valent (Kondo) limit, will be published elsewhere 
jnj . Using the Gutzwiller approximation, and minimiz- 
ing with respect to the mixing amplitudes, we find an 
expression for the ground state energy density £, which 
describes a renormalized d-band hybridized with effective 
/-levels: 
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where the k-sum extends over the Ut = Ud = Fermi sea 
(in a manner familiar from the Anderson lattice case [fio| , 
the Gutzwiller method respects Luttinger's theorem and 
leaves the Fermi volume unchanged) . The d-band kinetic 
energy renormalization factor qd is the same as in the 
Gutzwiller treatment |l3| of a Hubbard model with filling 
nd = Ti—nf /2. Since we assumed a dispersionless /-band, 
the g-factor for the /-electrons qf = (1 — n./)/(l — n.f/2) 
appears only in the renormalized hybridization amplitude 
v = ^Jqdqjv. The renormalized band structure contains 
the effective /-level e/; this shift is compensated by the 
last-but-one term of £ . The last term gives the d-d in- 
teraction energy: it contains Vd, the density of doubly 
occupied d-orbitals. £ still has to be optimized with re- 
spect tO Vd- 

Henceforth we assume a band with a constant den- 
sity of states p(e) — 1/W, lying in the interval e 6 
[— W/2, W/2]. Furthermore, we restrict ourselves to the 
weak hybridization regime v <C W , where a number of 
physical quantities shows an exponential Kondo-like be- 
haviour. One of them is the small deviation from integral 
/-valence. Minimizing with respect to Vd requires 
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which couples the non-integral valence to the Hubbard 
Ud- The solution is 
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In the expression for the optimized total energy density 
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the first line gives the energies of decoupled /-, and d- 
electrons. The coupling between the two subsystems is 
described by the last term which we identify as the Kondo 
energy of the Anderson-Hubbard model. It turns out to 
be proportional to 1 — n/. An analogous relationship 
between the valency and the Kondo temperature holds 
in the single impurity case Jl3] . 

In the above equations, n d = n — 1/2 is the v = 
value of the conduction electron density (per spin). q d is 
the g-factor taken with n d , and Ho{Ud) is the chemical 
potential of the Hubbard subsystem calculated (in the 
Gutzwiller approximation) for band filling n d . 

(^) and (|^) are formal relationships in the sense that 
they express the solution in terms of the optimized pa- 
rameters of the Hubbard subsystem, which have no sim- 
ple closed form for arbitrary Ud away from half- filling. 
We will give the results of small- Ud, and large- Ud, expan- 
sions. First, however, let us discuss the formal solution. 

(0) and (||) are the main results of the present paper. 
They say that switching on Ud influences the character- 
istic Kondo scale in two ways: through the prefactor and 
the exponent. The prefactor describes the correlation- 
induced narrowing of the d-band. The exponential factor 
can be written as exp (— W/ J c g) where we introduced the 
effective Kondo coupling 
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It differs from the Anderson lattice result [|I5|,|10| inas- 
much as /Zo (0) = {n d - l/2)W is replaced by the In- 
dependent chemical potential fJ.o(U d ). 

Our J e s differs from that of the corresponding impu- 
rity problem |l| in two respects. First, it has Av 2 in the 
numerator instead of 2v 2 : the "lattice enhancement of 
the Kondo effect" |l5| works also in the present model. 
Second, the denominator is also different, as we will dis- 
cuss shortly. 

In the weak coupling limit (Ud <§C W), the standard 
Hubbard model results 
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n d (l - n° d )(l - 2n° d )(Uj/W) into (@) and (§). The Kondo energy vanishes both at 

l d ( we run ou l °f conduction electrons) and at 
i° d — > 1/2 (because we are above the Brinkman-Rice 
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have to be replaced into (^) , or (j|) . We point it out that 
both in the prefactor and in the exponent, the appearance 
of U d acts to diminish the Kondo scale. 

It is worth to discuss the half-filled case n = 1 which 
by Luttinger's theorem arguments belongs to a correlated 
non-magnetic insulator. The second-order term in (|Io| ) 
vanishes at n d = 1/2. In fact, at half-filling no(Ud) = 
U d /2 holds for arbitrary U d . (Note that in the single- 
impurity result Q, Ud appears instead of our Ud/2). (|ll| ) 
reduces to the wellknown jl(| q d = 1 — (U d /2W) 2 which 
holds right up to U d = 2W. As U d —> 2W from below, 
the d-band undergoes a Brinkman-Rice transition, the 
system becomes integral valent (n/ = 1) and the Kondo 
effect is completely quenched. Within the Gutzwiller- 
type treatment of correlations, it is not unexpected that 
this should happen: the Kondo effect is due to the f—d 
hybridization which causes charge fluctuations in the d- 
shells. In the Brinkman-Rice scenario of the Mott tran- 
sition, these fluctuations get completely suppressed, and 
then so is the Kondo effect. — We should emphasize, 
though, that while the usual one-band Brinkman-Rice 
transition is a metal-insulator transition, here we are 
speaking about an insulator-insulator transition. For 
Ud < 2W, the effective d-band has a finite width and the 
renormalization of the hybridization matrix element gives 
rise to a modified Kondo gap. Calculating it requires a 
straightforward extension of the variational method to 
the n > 1 case. Having done this, the gap can be calcu- 
lated as the discontinuity of the chemical potential: 
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The result can be rewritten as A = 2qjv 2 /W , i.e., A 
can be interpreted as the renormalized hybridization gap. 
- For U d > 2W all polarity fluctuations cease and the 
system becomes a Mott insulator. 

Returning to the case of a general n d < 1/2, in the 
strong coupling limit (U d >• W) we replace 
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to the Ud — (Anderson lattice) case, we see that 
again, both the prefactor and the new term in the ex- 
ponent tend to reduce the Kondo scale. Taken together 
with the previous weak coupling result, wc conclude that 
the effective Kondo coupling is less than its Ud = 
value for all Ud > 0. For n slightly less than 1, the 
Anderson-Hubbard model predicts an even greater mass 
emhancement than the Anderson lattice, possibly ac- 
counting for findings like the heavy-fermion behaviour 
of Nd 2 _ 3; Ce 3; Cu04. 

It is curious that the small W/f/d-correction seen in 
(U|) changes sign as a function of the band filling n° d at 
rid = 1/3. If it were for the exponent alone, it might 
have been expected that at 1/3 < n d < 1/2, at suf- 
ficiently large U d , the Kondo scale increases towards its 
Ud — > oo value, and thus on the whole, its [/^-dependence 
in non-monotonic. However, our numerical results indi- 
cate that this tendency is counteracted by the stronger 
[/^-dependence of the prefactor and for reasonable pa- 
rameters, the Kondo scale is a monotonically decreasing 
function of Ud- We should mention, though, that this 
aspect of the behaviour is strongly influenced by the as- 
sumption about the form of the density of states. — All 
in all, our lattice result is qualitatively similar to the An- 
derson impurity results in the strong-coupling limit , 
while a discrepancy in the weak-coupling case H may be 
ascribable to the peculiarities of the Luttinger liquid. 

We should mention the limitations of our variational 
method. The use of the Gutzwillcr approximation is jus- 
tifiable only in the limit of infinite dimensionality but 
for many purposes, it should be a good approximation 
in three dimensions. In contrast, a number of previous 
studies dealt with a magnetic impurity imbedded in a 
one-dimensional system &H,||. An essential change of 
physics in going from one to two dimensions has been 
pointed out in [9|; three-dimensions should differ from 
both. - As for the limitations of the Ansatz (13), it 




shares with similar treatments of the Anderson 115 
and Kondo |l7| lattices the shortcoming of neglecting in- 
tersite correlations other than those following from the 
Fermi statistics of the projected Fermi sea. 

To summarize, we found that switching on the Hub- 
bard U d of the conduction electron sea tends to diminish 
the Kondo energy scale of the periodic Anderson model. 
This behaviour is just the opposite of that found for the 
Kondo-Hubbard lattice 101. We believe that this dis- 
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crepancy is understandable, and that the suppression of 
the Kondo scale in the Anderson-Hubbard lattice model 
arises from the fact that increasing Ud suppresses charge 
fluctuations, and therefore blocks the hybridization pro- 
cesses which would lead to the Kondo effect. Ud appears 
both in the exponent and in the prefactor. In the ex- 
ponent, we find the L^-dependent chemical potential of 
the d-band, while the prefactor reflects the correlation- 
induced narrowing of the conduction band, and it is 
sensitive to the nearness of the Brinkman- Rice transi- 
tion. Including the effects of a finite /-bandwidth, simi- 
lar conclusions should apply to the non-magnetic phases 
of such models as Tasaki's two-band model jL2j, or the 
Fu-Doniach model @]. It appears to us that the com- 
mon difficulty in many models of current interest is the 
description of Kondo compensation by a set of interact- 
ing electrons, and hence we hope that the relevance of 
the issues discussed here transcends the usual domain of 
/-electron physics. 
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